Crossover behaviors in liquid region of vortex states 
above a critical point caused by point defects 
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Vortex states in high-T c superconductors with point defects are studied by large-scale Monte 
Carlo simulations of the three-dimensional frustrated XY model. A critical point is observed on 
the first-order phase boundary between the vortex slush and vortex liquid phases. A step-like 
anomaly of the specific heat is detected in simulations of finite systems, similar to an experimental 
observation [F. Bouquet et at, Nature (London) 411, 448 (2001)]. However, it diminishes with 
increasing system size, and the number and size distribution of thermally-excited vortex loops show 
continuous behaviors around this anomaly. Therefore, the present study suggests a crossover rather 
than a thermodynamic phase transition above the critical point. 

PACS numbers: 74.25.Qt, 74.62.Dh, 74.25.Dw 



Intensive studies on vortex states of high-T c supercon- 
ductors have revealed complex structures of the phase 
diagram in the magnetic field along the c axis. Point 
defects such as oxygen deficiencies are essential for such 
complexity. The Bragg glass (BrG) 0, Q phase is sta- 
ble at low temperatures and low magnetic fields. This 
phase melts into the vortex liquid (VL) or vortex glass 
(VG) [3| phases as temperature T or magnetic field B 
increases, respectively. The BrG phase has a closed first- 
order boundary in the B-T phase diagram. 0, The 
vortex slush (VS) p phase locates between the VG and 
VL phases, and another first-order boundary between the 
VS and VL phases terminates at a critical point. jE0| 

The BrG-VL,J|! HElIll BrG-VG [illllEIS 
and VS-VL 0,0] phase boundaries have been observed 
experimentally. Similar phase diagrams have been ob- 
tained by numerical simulations including point defects, 
El El El though the stability of the VG phase 
is still controversial. p4l25| 

Recently Bouquet et al. [2(| observed a (5-function 
peak and a sharp step of the specific heat below and 
above a critical point, respectively. They argued that this 
step-like anomaly corresponds to a second-order phase 
transition related to a vortex- loop blowout. |27l l2q 

In the present article, we investigate the phase diagram 
in the vicinity of the critical point of the VS-VL phase 
boundary. As displayed in Fig.^ we find that the critical 
point is located at e « 0.15, which is just above the 
maximum value of e in our previous work. 23] A step-like 
anomaly of the specific heat is observed above the critical 
point in finite systems. However, the size dependence of 
this anomaly reveals that this step does not indicate a 
thermodynamic phase transition but a crossover. 

In the present study we use the same model as investi- 
gated in our previous article on vortex states with point 
defects, |23l namely the three-dimensional anisotropic, 
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FIG. 1: Phase diagram of the model Q in the plane of the 
strength of point defects (e) and temperature based on a finite 
system with L a — Lb = 50 and L c = 40. Part of this phase 
diagram (for e < 0.14) has already been published. (23] The 
first-order VS-VL phase boundary terminates at the critical 
point (full circle), and the step of the specific heat (dotted 
line) is observed above the critical point. The VG-VS phase 
boundary for e > 0.14 is out of the scope in the present study. 



frustrated XY model on a simple cubic lattice, [2(1 Isfl ] 



H 



Am — 



2tt 



i - 4>j Ajj ) 

m $n ) ) 

- I" A^-drM, B = VxA, 

Ji 



i.j£ab plane 

-- y 



in, -n \\c axis 



(i) 

(2) 



*Also at Lyman Lab. of Physics, Harvard University, Cambridge, 
MA 02138; Electronic address: nonomura.yoshihiko@nims.go.jp 



with the periodic boundary condition in all directions. A 
uniform magnetic field B is applied along the c axis, and 
the averaged number of flux lines per plaquette is /. A 
point defect is introduced as the plaquette surrounded by 
four weaker couplings and distributed randomly in the ab 
plane with probability p. Couplings are Jy = (1 — e) J 
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FIG. 2: Temperature dependence of (a) the specific heat, (b) 
the density of dislocations in the ab plane and (c) the ratio of 
point defects which trap flux lines for e = 0.14, 0.15, 0.16, 0.18 
and 0.2 (from bottom to top). The origin of the quantities 
are varied for each e in (a) and (b). 



(0 < e < 1) on the point defects, and Jy = J elsewhere. 
As in our previous article, we fix the parameters T = 20, 
/ = 1/25 and p = 0.003, and vary the strength of point 
defects e and temperature T. Configurations of defects 
are also the same as in our previous work. Monte Carlo 
simulations are started from a very high temperature, 
and systems are gradually cooled down. Typical Monte 
Carlo steps are 4 ~ 14 x 10 7 for equilibration, and 2 ~ 
16 x 10 7 for measurement. 



First, we concentrate on the system with L a — Lb = 50 
and L c — 40. Temperature dependence of the specific 
heat C for 0.14 < e < 0.2 is displayed in Fig.|2». At 
e = 0.14, C has a (^-function peak at T s \ = 0.0760 J/kB- 
This peak becomes very small at e = 0.15, and is replaced 
by a step for e > 0.16. This result suggests that the 
critical point locates at e » 0.15. The phase diagram 
obtained from this finite system is summarized in Fig. ^ 

Consistent behaviors are seen in other quantities. Tem- 
perature dependence of the density of dislocations in the 
ab plane pd and the ratio of point defects which trap flux 
lines 7"t r are given in Figs.JSJb) and[5Jc), respectively. The 
quantity p^ characterizes 0, |23| the VS-VL phase transi- 
tion. At e — 0.14, pa jumps between T = 0.072 J / k& and 
0.076J//cb, where C sharply increases. This jump be- 
comes much smaller at e = 0.15 and invisible for e > 0.16. 
Keeping possible statistical errors in mind, a tiny hump 
and dip seem to exist at the temperature region T* where 
the specific heat has a step-like anomaly for e > 0.15, but 
the hump and dip are one order smaller than the jump of 
Pd related to the VS-VL phase transition, and therefore 
considered to be finite-size effects. 

At e = 0.14, the quantity r tr gradually increases from 
T = 0.070J//cb as approaching the transition tempera- 
ture T s i, because the short-range order of flux lines in 
the VS phase becomes fragile as temperature increases. 
Then, it rapidly decreases just below T s \, because the 
mobility of flux lines becomes larger in the coexisting re- 
gion of the first-order phase transition. This anomaly 
becomes obscure at e = 0.15 and disappears for e > 0.16. 
Although a tiny hump and dip also seem to exist at T* 
for e > 0.15, they have an opposite behavior to the cusp 
at e — 0.14; r tr increases at T* and drops just above 
T*, and therefore have a different origin from the cusp 
related to the VS-VL phase transition. 

Next, in order to explore whether this step-like 
anomaly of the specific heat is a thermodynamic phase 
transition or not, its size dependence is investigated. 
Since flux lines in B || c are one-dimensional objects 
along the c axis, physical quantities mainly depend on 
L c . |3l| Even in the vortex liquid region in pure systems, 
the displacement of flux lines along the ab plane is smaller 

1/2 

than the random-walk value ~ L c due to the blocking 
by other flux lines. In the present case, this displacement 
would be even smaller due to the pinning by defects. We 
thus fix L a = Lb = 50 and change L c . Temperature de- 
pendence of the specific heat for L c = 40, 56 and 80 at 
e = 0.2 is displayed in Fig. Efa). As L c increases, the 
height of the step of the specific heat becomes smaller. 
Since it is clearly proportional to 1/L C (Fig. Efb)), this 
step is expected to vanish in the thermodynamic limit. 
The anomaly of the specific heat above the critical point 
of the VS-VL phase boundary is therefore a crossover 
rather than a thermodynamic phase transition. 

Finally, the number and size distribution of thermally- 
excited vortex loops are investigated in order to check 
whether this anomaly of the specific heat is related to 
the vortex- loop blowout [271 |28[ or not. Temperature 
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FIG. 3: (a) Temperature dependence of the specific heat at 
e = 0.2 for L c — 40, 56 and 80. Its origin is varied for each 
e. (b) c-axis size dependence of the height of the step in (a). 
The straight line is drawn by the least-squares fitting. 



dependence of the number of vortex loops for 0.14 < e < 
0.2 is displayed in Fig. 0Ja). There exist no vortex loops 
including pancake vortices at these parameter regions. 
Size distributions of Josephson loops for e = 0.14 and 0.2 
are shown in Figs.^b) and^Jc), respectively, ^ij^ 

Although the number of vortex loops monotonically 
increases as temperature increases, no symptom of the 
vortex-loop blowout is seen around T s \ or T* for each e. 
This is also the case for the size distribution of Josephson 
loops. In Figs.^b) and^c), the number of Josephson 
loops Nj(Lj) with a given perimeter Lj is plotted versus 
Lj in a semi-log scale, and the exponential size distri- 
bution of Nj(Lj) is clearly observed. There exists no 
qualitative difference in the behaviors below and above 
the critical point, and the size distribution changes con- 
tinuously as temperature is varied around T s \ or T*. 

The vortex-loop blowout is characterized by the abrupt 
increase of the number of vortex loops and the change 
of the size distribution from an exponential decay to a 
power-law one at the transition temperature, and the 
proliferation of vortex loops as large as the system size 
above the transition temperature. |27ll28| None of these 
behaviors are observed around T* in the present simu- 
lations, which casts serious doubt on the argument that 
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FIG. 4: (a) Temperature dependence of the number of vortex 
loops for e = 0.14, 0.15, 0.16, 0.18 and 0.2 (from bottom 
to top). Its origin is varied for each e. Size distributions of 
Josephson loops for (b) e = 0.14 and (c) e = 0.2 around the 
anomalies of the specific heat are drawn in a semi-log scale. 



the step of the specific heat is related to the vortex-loop 
blowout. [26( In this argument, the vortex-loop blowout 
scenario was borrowed in order to explain the observa- 
tion that the value of the specific heat above the step-like 
anomaly C+ is larger than that below the anomaly C_. 
However, this property is quite general in vortex states, 
and actually occurs even at the melting transition in pure 
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systems. |3(j In vortex states, the density of vortices 
increases as temperature increases owing to excitation 
of vortex loops and wandering of flux lines, and conse- 
quently the relation C+ > C_ holds. The vortex-loop 
blowout is merely a special case of this general property. 

It remains unsolved why a step of the specific heat is 
clearly observed in experiments even though the sample 
size is much larger than the system size in our simula- 
tions. Several possibilities are: (a) The height of the step 
of C might saturate at larger L c , and our maximum sys- 
tem size L c = 80 would be still small, (b) Our model 
JTJ might be still too simple to describe the step of C, 
e.g. the inhomogeneity of the strength of point defects 
is not taken into account, (c) The step of C observed 
experimentally might be a non-equilibrium effect. 

In summary, vortex states of high-T c superconductors 
with point defects in B || c has been investigated by 
large-scale Monte Carlo simulations. The critical point of 
the first-order phase boundary between the vortex slush 
and vortex liquid phases is identified as the vanishing 
point of the S- function peak of the specific heat. The 



shape of the step-like anomaly of the specific heat above 
the critical point is similar to that of the experiment by 
Bouquet et al. As the system size along the c axis L c 
increases, the height of this step seems proportional to 
1/L C , which suggests that this anomaly is a crossover 
rather than a thermodynamic phase transition. Consis- 
tent behaviors are observed in the temperature depen- 
dence of the density of dislocations in the ab plane and 
the ratio of point defects which trap flux lines. There 
is no discontinuity in the number and size distribution 
of vortex loops around the step-like anomaly of the spe- 
cific heat, which clearly indicates that this anomaly is 
not related to the vortex- loop blowout. 
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Sports, Science and Technology, Japan. This study is 
partially supported by the Priority Grant No. 14038240 
from the same Ministry. 



[1] T. Nattermann, Phys. Rev. Lett. 64, 2454 (1990). 

[2] T. Giamarchi and P. Le Doussal, Phys. Rev. Lett. 72, 

1530 (1994); Phys. Rev. B 52, 1242 (1995). 
[3] M. P. A. Fisher, Phys. Rev. Lett. 62, 1415 (1989); 

D. S. Fisher, M. P. A. Fisher, and D. A. Huse, Phys. 

Rev. B 43, 130 (1991). 
[4] D. Ertas. and D. R. Nelson, Physica (Amsterdam) 272C, 

79 (1996). 

[5] T. Giamarchi and P. Le Doussal, Phys. Rev. B 55, 6577 
(1997). 

[6] T. K. Worthington et al, Phys. Rev. B 46, 11854 (1992). 
[7] J. Kierfeld and V. Vinokur, Phys. Rev. B 61, R14928 
(2000). 

[8] H. Safar et al, Phys. Rev. Lett. 69, 824 (1992). 

[9] W. K. Kwok et al, Phys. Rev. Lett. 72, 1092 (1994). 
[10] E. Zeldov et al, Nature (London) 375, 373 (1995). 
[11] U. Welp et al., Phys. Rev. Lett. 76, 4809 (1996). 
[12] A. Schilling et al, Nature (London) 382, 791 (1996). 
[13] R. Cubitt et al., Nature (London) 365, 407 (1993). 
[14] H. Safar et al., Phys. Rev. B 52, 6211 (1995). 
[15] B. Khaykovich et al., Phys. Rev. B 56, R517 (1997). 
[16] T. Nishizaki et al, Phys. Rev. B 58, 11 169 (1998). 
[17] T. Nishizaki, K. Shibata, and N. Kobayashi, Physica 
(Amsterdam) 341C-348C, 957 (2000); K. Shibata et al, 
Phys. Rev. B 66, 214518 (2002). 
[18] H. H. Wen et al., Phys. Rev. B 64, 054507 (2001). 
[19] M. J. P. Gingras and D. A. Huse, Phys. Rev. B 53, 15193 
(1996). 



[20] S. Ryu et al., Phys. Rev. Lett. 77, 2300 (1996). 

[21] A. van Otterlo, R. T. Scalettar, and G. T. Zimanyi, Phys. 

Rev. Lett. 81, 1497 (1998). 
[22] R. Sugano et al., Physica (Amsterdam) 284B-288B, 803 

(2000) , 341C-348C, 1113 (2000). 

[23] Y. Nonomura and X. Hu, Phys. Rev. Lett. 86, 5140 

(2001) . 

[24] P. Olsson and S. Teitel, Phys. Rev. Lett. 87, 137001 
(2001). 

[25] A. Vestergren, J. Lidmar, and M. Wallin, Phys. Rev. 

Lett. 88, 117004 (2002). 
[26] F. Bouquet et al, Nature (London) 411, 448 (2001). 
[27] Z. Tesanovic, Phys. Rev. B 51, 16204 (1995), 59, 6449 

(1999). 

[28] A. K. Nguyen and A. Sudb0, Europhys. Lett. 46, 780 

(1999); Phys. Rev. B 60, 15307 (1999). 
[29] Y.-H. Li and S. Teitel, Phys. Rev. Lett. 66, 3301 (1991); 

Phys. Rev. B 47, 359 (1993). 
[30] X. Hu, S. Miyashita, and M. Tachiki, Phys. Rev. Lett. 

79, 3498 (1997); Phys. Rev. B 58, 3438 (1998). 
[31] Y. Nonomura, X. Hu, and M. Tachiki, Phys. Rev. B 59, 

R11657 (1999). 

[32] The effect of intersections of vortex loops is negligible 
in the present parameter regions, because the averaged 
number of vortex loops in each ab plane is less than unity, 
in contrast to the case of Ref . |33(l . 

[33] P. Olsson, Europhys. Lett. 58, 705 (2002). 



